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Abstract 

The contribution of gravitational neutrino oscillations to the solar neutrino 
problem is studied by constructing the Dirac Hamiltonian and calculating 
the corresponding dynamical phase in the vicinity of the Sun in a non- 
Riemann background Kerr space-time with torsion and non-metricity. We 
show that certain components of non-metricity and the axial as well as 
non-axial components of torsion may contribute to neutrino oscillations. 
We also note that the rotation of the Sun may cause a suppression of 
transitions among neutrinos. However, the observed solar neutrino deficit 
could not be explained by any of these effects because they are of the order 
of Planck scales. 
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1 Introduction 



Neutrinos always attracted a lot of attention in high energy physics pQ. A major 
problem of interest at present is the solar neutrino problem. The Sun is a strong 
source of electron neutrinos v e because of the thermonuclear reactions taking 
place in its core. According to the standard solar model, the number of v e to be 
emitted from the Sun can be predicted. At the same time, the flux of electron 
neutrinos coming from the Sun can be measured on earth. The measured amount 
of u e is approximately one third of the predicted amount. Essentially, this is 
the so-called solar neutrino problem. One well known solution to this problem 
is provided by the assumption of neutrino oscillations PQ,[2]- Briefly stated, the 
neutrino oscillations imply that the electron neutrinos coming out of the Sun may 
be converted to other neutrino species, muon and tau v T during their journey 
towards the earth, assuming neutrinos to have a mass whereas the standard 
electroweak model asserts zero mass for them. It should also be noted that all 
the above arguments have been cast in Minkowski space-time. However, we know 
that we live in a curved space-time - perhaps even in a curved space-time with 
torsion and non-metricity. Therefore, in more recent years, physicists have turned 
their attention to specifically gravitational contributions to neutrino oscillations - 
see [3] , |H > , [H] > [Zj and references therein. We recently investigated the effects of 
space-time torsion on neutrino oscillations |8 r and see also 0,(101 • The essence 
of this work is to calculate the dynamical phase of neutrinos, by finding the form 
of the Hamiltonian, H, from the Dirac equation in a non-Riemannian space-time. 
The phase then follows from the formula 



where tp is a Dirac 4-spinor and H is a 4 x 4 matrix. The Hamiltonian H will 
depend, for example, on momentum p, and this is expressed not as a differential 
operator but simply as a vector. 1 In this note we investigate within the same 
approach the possible effects of space-time non-metricity on neutrino oscillations. 



Space-time is denoted by the triple {M, g, V} where M is a 4-dimensional dif- 
ferentiable manifold, equipped with a Lorentzian metric g which is a (0,2)-type 
covariant, symmetric, non-degenerate tensor and V is a connection which de- 
fines parallel transport of vectors (or more generally tensors). We shall give a 
coordinate system set up at a point, p £ M, by coordinate functions (or inde- 
pendent variables) {x a (p)}, a = 0, 1,2,3. This coordinate system forms a set of 
natural (or coordinate) reference frame at p as {-J^s(p)}, with shorthand notation 
d a = q^jj. This natural reference frame is a basis vector set for the tangent space 

1 The exponential of e _7r J Hdt is denned by its power series expansion. 





2 Space-Time Geometry 



l 



dx a {^) = ia,dx a = S% . (2) 



at p, denoted by T P (M). Similarly, differentials {dx a (p)} of coordinate functions 
{x a (p)} at p, form a natural (or coordinate) reference co- frame in the co-tangent 
space at p, denoted by T*(M). Interior product of the basis vectors with the 
basis co-vectors is defined by the Kroenecker symbol: 

d_ 

In general, any set of linearly independent vectors in tangent space, T p (M), can be 
taken as basis vectors and these vectors can be ort honor malized by, for example, 
the Gram-Schmidt process. We denote a set like this by {X a }, a = 0, 1, 2, 3 and 
call it an orthonormal reference frame. In this case the metric defined on M 
satisfies the relation 

g(X a ,X b ) = rj ab (3) 

where 7] ab is known as the Minkowski metric which is a matrix whose diagonal 
terms are -1,1,1,1 and off-diagonal terms are zero. The basis set dual to the 
orthonormal reference frame are denoted by {e a }, a = 0,1,2,3, and called the 
orthonormal reference co-frame. {X a } and its dual {e a } satisfy the following set 
of equalities that is another manifestation of eqn.Q: 

e a (X b )=z Xb (e a ) = 5 a b . (4) 



Here we adhere to the following conventions: indices denoted by Greek letters a, 
/3, • • • = 0, 1, 2, 3 and /1, u, • • • = 1, 2, 3 are holonomic or coordinate indices, a, 
b, • • • = 0, 1, 2, 3 and i, j, ■ ■ ■ = 1, 2, 3 are anholonomic or frame indices. In 
terms of the local coordinate frame d a (p), the orthonormal frame X a (p) can be 
expanded via the so-called vierbein (or tetrad) h a a (p) as 

X a (p) = h a a (p)d a (p) . (5) 

In order for X a to serve as an anholonomic basis, the h a a (p) are required to be 
non-degenerate, i.e., deth a a (p) 7^ 0. In T*(M) an orthonormal co- frame e a (p) can 
be expanded in terms of the local coordinate co- frame dx a (p) as 

e\p) = h%(p)dx?(p) . (6) 

The inverse vierbein h b p{jp) have to be non-degenerate as well. Moreover, the 
duality of the frame and the co-frame requires for the vierbein and its inverse to 
satisfy 

ix a e b = h\{p)h b a {p) = 6 b a . (7) 

We set the space-time orientation by the choice €0123 = 1- The non-metricity 1- 
forms, torsion 2-forms and curvature 2-forms are defined by the Cartan structure 
equations 

2Qab = -D Vab := A ab + A ba , (8) 
T a = De a := de a + A\ A e b , (9) 
R\ = DA a b := dA a b + A a c A A c b . (10) 
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d, D, i a , * denote the exterior derivative, the covariant exterior derivative, 
the interior derivative and the Hodge star operator, respectively. The linear 
connection 1-forms can be decomposed in a unique way according to [TT] : 

K\ = uj a b + K a b + q\ + Q\ (11) 

where uo a b are the Levi-Civita connection 1-forms 

uj a b A e b = -de a , (12) 

K a b are the contortion 1-forms 

K\ A e b = T a (13) 

and q a b are the anti-symmetric tensor 1-forms 

q a b = -(i a Qbc)e c + Mac)e c . (14) 

It is cumbersome to take into account all components of non-metricity and 
torsion in gravitational models. Therefore we will be content with dealing only 
with certain irreducible parts of them to gain physical insight. The irreducible 
decompositions of torsion and non-metricity invariant under the Lorentz group 
are summarily given below. For details one may consult Ref. [T2] . The non- 
metricity 1-forms Q ab can be split into their trace-free Q ab and the trace parts 
as 

Qab = Qab + \vabQ (15) 

where the Weyl 1-form Q = Q a a and rj ab Q ab = 0. Let us define 

A& := i a Q a b , A:=A a e a , 

:= *(Q ab ^e a ), Q:=e b AQ b , Q a :=Q a -^ a Q (16) 

as to use them in the decomposition of Q ab as 

Qab = Q^b + Qa\ + Qab + Qab ( 17 ) 



where 



Qi? = ^*(e a AQ b + e b AQ a ) (18) 

Qab = l(A a e b + A b e a -^ Vab A) (19) 

Q { £ = \vabQ (20) 

QS = Qab-Q { S-Q { 2-Q { ab ] ■ (21) 
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We have i a Q§ = i a Q^ = 0, V ab Q% = V " b Q® = jfbQ® = 0, and 
e a A Qoj} = 0. In a similar way the irreducible decomposition of T a 's invariant 
under the Lorentz group are given in terms of 



so that 



where 



a = i a T\ a = e a AT a (22) 



rpa _ rpa(l) _|_ rpa(2) _|_ rpa(3) {23) 



T a{2) = ^Aa, (24) 

T a{3) = -i a a , (25) 
3 

t a ■= r p a (X) — r £ a r £ a ^) rjia(3) (26) 

Here i a t a = t a T a{3) =0, e a A t a = e a A T a{2) = 0. To give the contortion 
components in terms of the irreducible components of torsion, we firstly write 

2K ab = i a T b - i b T a - (i a i b T c )e c (27) 

from (fT3~|) and then substituting (J2*3~]) into above we find 

2K ab = i a t b - i b t a - (i a ibtc)e c 

+ g(e a A i b a - e b A i a a) 

2 1 
+-(v&cr) - -(z a z^ c a)e c . (28) 

In components if a6 = if Ciafe e c , t a = \t bC)CL e bc , a = F a e a , cr = ^a abc e abc this 
becomes 

K c ab — (tac,b tbc,a tab,c) 

+ \(F b r] ac - F a r] bc ) - \a abc . (29) 

3 O 

3 Hamiltonian of a Dirac particle in arbitrary 
space-times 

The Dirac equation in a non-Riemannian space-time with torsion and non-metricity 
is written as [TH] . [H] . [TH] 

*7A^ + M*1^ = (30) 
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in terms of the Clifford algebra Clz,i~ valued 1-forms 7 = 7 a e a and M = We 
use the following Dirac matrices 



.10] t . ( (T ?; 







where a 1 are the Pauli matrices, ip is a 4-component complex valued Dirac spinor 
whose covariant exterior derivative is given explicitly by 

DiP = dijj + ^A [afe W + \Qip (31) 



where 



^afc = ^[7a,7b] (32) 



are the spin generators of the Lorentz group. We write it out explicitly as 

Dip = dt^ + d^dx^ + ^n^a^ + ^Q^ 

= h 6 c e c d t ?P + h\e c d^ + l -e c Sl CAh o ab i> + ^e c Q^ (33) 

where Auw := fi a t = Q Cta i,e c anti-symmetric part of the full connection 1-form 
and Q = Q a e a and using *7 = r y a *e a and the identity *e a A e b = —r) ab *l calculate 

*7 A Dip = (-hK^dtiP - W cl c d^ - ^O c , o6T V V - \Qd c ^j *1 • (34) 
Putting this into (J3Uj) we obtain 



h° cl c d t ip = -W cl c d^ + Mip- -n cMl c a ab tP - - A Q cl c tp. (35) 



We multiply this from left by 

^hKry^^frW) (36) 

where 

b 2 :=(h\) 2 + h\h^ . (37) 
When we compare the result with the Schrodinger equation 

ih^ = Hip, (38) 

we deduce the Dirac Hamiltonian matrix |Zl (HI UHl UEl EIZI EIHl EIHl EHj 

h = ^hKh^m.-^hK^ 

+ ^h° d {l c , abl d l c a ab + |^ 6 aQ 6 7 a 7 6 - (39) 



The right hand side of ()39|) need not be a hermitian matrix in general; e.g. if 
h°i ^ 0, then the mass term contains an anti-hermitian part such as 

H = H + iHt (40) 

where Hq = Ho and Hi = H\. However, the decomposition (|4()jl is frame 
dependent. That is we can always find a local Lorentz frame in which Hamiltonian 
is fully hermitian fTf\ ITS] . First we can get rid of the anti-hermitian part of the 
mass term by diagonalizing the matrix h a a via a frame transformation 

d a (x) -> df3(x)L~ ll3 a (x) 

g a p(x) -> g^ 5 (x)L~ 1 ' y a L' 1 %. (41) 

Thus 2 

h a a (x) -> f(x)8° (42) 

where x stands for x a and f(x) is composed of h a a (x). Under this change (}3*9"j) 
goes over to 

H -> H = cf 1 (x)-f°-f i ihd i -imc 2 f 2 (x) 1 ° 

+ihch{x)^ abl Q 1 c a ab + inc/ 4 (x)Q 67 V (43) 

where are composed of h a a (x). Putting in the definition 

^c,afe = £abcd.S d (44) 

and using the identity 

7 v 6c = I^y _ I^y - I e ^ 7d75 (45) 

where 75 = 70717273, the Hamiltonian matrix becomes 

H = cMxh^thd- - tmc 2 f 2 (xh° + ihcNaix)^ 11 + ihcf B (x)S a >y j a >y s (46) 
where we introduced 

N a :=f 3 (x)Q\ a + U(x)Q a . 
If we now define the canonical momenta 

Pi ■= -ih ( ft + ^rr) (47) 



and assume 



hi*) 



Pi =pi, (4£ 



2 L e SO+(l,3) where SO+(l,3) is special orthochronous Lorentz group. 
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PUj) takes the form 

H = /i(x)cp J 7o7i + imc 2 f 2 (x)'j + ihcf 5 (x)S a 'j 'j a 'j 5 - ihcN (x) . (49) 

In order eliminate the last term in (jH^ one may further perform a locally unitary 
transformation 

ip(x) -»■ U + (x)ip(x) , H -> U + (x)HU(x) (50) 

and obtain 

H ^ H = f\(x)cp l U + (x) loli U(x) + imc 2 f 2 (x)U + (x)-y U(x) 
+thcf 5 (x)S a U + (x) l0lal5 U(x) 

-ihcih^U+ix^d'Uix) + JVo(z)] • (51) 
Under the following solvable matrix equation 

U+Wmi&Uix) = (52) 

we give the final form of our hermitian Hamiltonian matrix (up to a sign) by the 
expression 

H = fi{x)cp l j ji + imc 2 / 2 (x)7q + ihcf 5 (x)S a j j a j 5 . (53) 

4 Neutrino oscillations in the Kerr background 

Here we construct the Hamiltonian matrix of a Dirac particle (i.e. a massive 
neutrino) of mass m in the background space-time geometry of a heavy, slowly 
rotating body of mass M such as the Sun. Its exterior gravitational field will be 
described by weak constant, uniform torsion and non-metricity fields, together 
with the Kerr metric |21j : 

/ 2MGr\ p 2 

ds 2 = - 1 5-^- cdt®cdt+ ^-dr <g> dr + P 2 d9 <g> d6> 

y cr^cr y A 

/ 2 a 2 2MGa 2 r . , \ . 2 „ , , 4MGar . 2 „ , 

+ ( r + H ^ — sin 2 8 sm 2 6 dip ® dip sm 2 6 dt ® dip (54) 

V <r c 4 p y ' ' c z p 2 

where A = r 2 - + (f ) 2 , p 2 = r 2 + (f ) 2 cos 2 6 ,a = ^ = \R 2 u. The Sun is 
assumed a uniform sphere of radius R. M, J and u are the mass, angular momen- 
tum and angular velocity of the Sun, respectively. We choose the orthonormal 
co-frame 



e 2 = pd0 , e z = ^tt r 2 +(^ 2 )dip-adt\ (55) 
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and using the definitions 

111 

de a + u a b A e b = <^> to ab = --i a de b + -i b de a + -(2 a 2 6 cfe c )e c , (56) 

calculate the Levi-Civita connection 1-forms 

MG[r 2 - (|) 2 cos 2 9] [(Mf - r)p 2 - ^f^]a sin 2 9 

uj 1 - -f dt H - — - dip, 

p 4 c p 4 c 

2 2MGra cos 9 , A(^) 2 sin 2 9 - (r 2 + (^) 2 ) 2 

u>\ = — dt+^^ j± \cL±_ cos 0dtp, 

p 4 c 2 p 4 

\fKa sin 9 cos 

w 2 = j dp, 

p A c 



i v Ar sin 
3 = - dip, 

p 2 



VKa cos 9 ar sin u , 

uj , = d9 = — ar, 

P 2 c ^Ap 2 c 



, a 2 sin 6> cos 6> , r\/A ,„ . , 

= = dr--^—d9. 57 

2 p 2 ^ 2 P 2 V ; 

To simplify the discussions, we consider only the motion of massive neutrinos 
restricted to the equatorial plane of the Sun. Thus we set 9 = ir/2 and d9 = 
0. Furthermore, since the Sun rotates very slowly [uo ~ 3 x 1CT 6 (rad/s)} we 
approximate the metric functions. Therefore, in reasonably far away distances 
from the Sun, the restricted line element will be taken as 

, 2 / 2MG\ , , / 2MGV 1 , 

ds 2 ~ - 1 — )cdt®cdt+ [1 - — dr®dr 

V err ) \ c 2 r J 

a MG 

+r dip ® r dip — 4 — — - cdt £g> r dp . (58) 
c c ' r 

We also write the orthonormal co- frame approximately up to O (^) as 

e ° = f c dt - —dp , e 1 = 4^ , e 2 = , e 3 = --di + rdp (59) 
c / r 

where 

p^ err 
The inverses of these relations to the same order of approximation are 

cdt = - f e° + — e 3 , rfr = /e 1 , dO = , dy? = -^- e ° + -e 3 (61) 
/ rc ' jr 2 c r 



which give 



h\ = - r /> 6 3 = -, h\ = f, h\ = -^- 2 , h\ = - (62) 



with all the other components neglected. To this order of approximation (|57|1 
gives 

^fe° + ^e 3 , u; 03 ^— 2 e\ u 3l ~ + L 3 (63) 

with the remaining ones neglected. Then the Hamiltonian matrix (|39|) reads 

# ~ / 2 q» r 7o7i + /cp«p7o73 + ifmc 2 ^ ~ -=hcf f'loli 
3 

+-i/zc/^ a 7o7 5 7a - ihcfN a ^ a 
af 3 iahpf iamcf 2 

Pr737l H o 7371 73 

t At t 

iakf 3iahf 2 ictfif 2 

7o7275 H o 5 737a75 H iV 7 3 7a (64) 



2r 2r r 



where 
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Pr 


:= -zft— + ~ 
ar r 


Pe 


:= 




z?i d 


P<p 


r dip 


N a 


'■= -^{t b a,b + F a 
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(65) 
(66) 
(67) 

(68) 



Note that the contributions of axial components of torsion are given by S a while 
certain components of non-metricity and the non-axial components of torsion 
occur only in N a and the rotation effects are given in terms of the parameter a. 
We rewrite the Hamiltonian 4x4 matrix in terms of 2 x 2 matrices as follows: 



21 -'-'22 



with 



af 3 

H n = fmc 2 + iA + Bo x + (C - i— p r )a 2 + Da 3 



r 

af 3 

#22 = -fmc 2 + iA + Ba x + (C - i—p r )^2 + Da 3 

r 

CL I TI C 

Hu = F + (f 2 cp r + iG)a 1 +iHa 2 + (fcp v + — + iK)a 3 

r 

H 2 i = F+(f 2 cp r + tG)a 1 + tHa 2 + (fcp ip J — + iK)a 3 (70) 



where we set 



A ~ -hcfN + ^N 3 
r 
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_, 3. _, ahf 2 „ T 

2 r 

3 ofi f 2 

C ~ -^-^-(l + r/' + rJVi) 

3 3ahf 2 
D ~ -hcS 3 --^-S 

3 3aft/ 2 
2 ° 2r 3 
r W 3aft/ 2 

2r 

X ~ -^c/iVs + ^iV . (71) 



The way we approach the solar neutrino problem starts by writing down the 
Dirac equation in a rotating, axially symmetric background space-time geometry 
and finding phases corresponding to neutrino mass eigenstates, then finally calcu- 
lating the phase differences among them. There are two cases of special interest: 
the azimuthal motion and the radial motion. The analysis of the azimuthal mo- 
tion with p = (p r ,Pe,Pip) = (0, 0,p) yields for ultrarelativistic neutrinos, for which 
pc ~ E and cdt ~ Rdip, the phase for the spin up state 

^(/* + ^ + ^+,, + *))^ (72) 

and similarly for the phase of the spin down state 
where 

A v ~ B 2 + C 2 + D 2 + F 2 + G 2 + H 2 + 2(DF + BH - CG) . (74) 

These phases alone do not have an absolute meaning; the quantities relevant for 
the interference pattern at the observation point of the neutrinos are the phase 
differences A<3> = $ 2 — $x where $i and $2 are the absolute phases of the neutrino 
mass eigenstates v\ and u 2 . It is thus seen from equations (J72j) and (|73*|l that the 
phase differences can have explicit dependence on non-metricity in the case of 
opposite spin polarizations of mass eigenstates for the azimuthal motion via (|71j) : 

a *^-*-(Htt 2 ^^ (75) 
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where Am 2 = — mi 2 . 

The Hamiltonian for the radial motion on the other hand is obtained by the 
assumption p — (p, 0,0). In this case with the further assumptions pc ~ E and 
cdt ~ dr, the phases appropriate to the spin up and spin down particles are, 
respectively, 

*' = ±j(fE + ^ + yfc + i{A + G)y, (77) 
& = L[(f E + 1 ^-jA~ r + i(A + G))dr (76 

where 



hcJ \ 2E v v 7 



a , ^ r n9 , ,-, ^ amH amG , 
A r ~ {D-Hf + (B + F + 2 + C + iT 

a 2 / 4 , , 2 2m/ 2 amG, 



(mc-/j9) 2 + — ^(mc-/p)(C + K ). (79) 

In this case the relevant phase differences depending on non-metricity via N a and 
rotation via a come from the opposite spin polarization states 

/\ ID, ^ C*~* 2 f I 

A * = *5-*! = W Ar + ^/v^*' < 81 > 

We point out that A r = ReA r + 2lmA r implies \/A^ — ot + i(3 and hence the 
rotation of the Sun would suppress the transitions among the neutrinos via the 
phase difference equations (f£>0 J) .(|%T |) in opposite spin polarizations. 



5 Conclusion 

We have here extended our recent study of gravitationally induced neutrino oscil- 
lations jH] by including the effects of rotation of the Sun, space-time non-metricity 
and as well as components of torsion other than the axial ones. The rotation of 
the Sun implies a damping of neutrino oscillations. However, this result is frame 
dependent as we explained in Sect 121 in general. We have shown that there are 
contributions coming from non-axial components of spacetime torsion and def- 
inite components of spacetime non-metricity depending on the polarizations of 
the spin states of the mass eigenstates. If we set the rotation parameter a = 0, 
then (fTH|) gives 

^ = he ((So + SO 2 + (S a - 2 -fN,f + (S 3 + |/JV 2 ) 2 ) (82) 
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which means that there is no suppression among the neutrinos and only N2 and 
N 3 components of N a contribute to the oscillations. If we further set N a = 0, 
we reach agreement with our previous results in 0. It should be clear that the 
above scheme only works if the neutrino masses are different from in each other 
and hence, in general, different from zero. This means there are right-handed 
neutrinos as well as left-handed ones which, however, must interact with matter 
very weakly as they have not yet been observed. Finally, we note that all the 
possible contributions discussed here so far would be of the order of Planck scales, 
and hence do not suffice to account for the observed solar neutrino deficit. 
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